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Abstract 

We present a measurement of inclusive J/xj/ production in p-Pb collisions at = 5.02 TeV as a 
function of the centrality of the collision, as estimated from the energy deposited in the Zero Degree 
Calorimeters. The measurement is performed with the ALICE detector down to zero transverse 
momentum, pj, in the backward (—4.46 < jcms < —2.96) and forward (2.03 < jcms < 3.53) rapidity 
intervals in the dimuon decay channel and in the mid-rapidity region (—1.37 < jcms < 0.43) in the 
dielectron decay channel. The backward and forward rapidity intervals correspond to the Pb-going 
and p-going direction, respectively. The px-differential J/ 1 // production cross section at backward and 
forward rapidity is measured for several centrality classes, together with the corresponding average 
Pj and pj values. The nuclear modification factor is presented as a function of centrality for the 
three rapidity intervals, and as a function of pj for several centrality classes at backward and forward 
rapidity. At mid- and forward rapidity, the J/\j/ yield is suppressed up to 40% compared to that 
in pp interactions scaled by the number of binary collisions. The degree of suppression increases 
towards central p-Pb collisions at forward rapidity, and with decreasing pj of the J/xj/. At backward 
rapidity, the nuclear modification factor is compatible with unity within the total uncertainties, with 
an increasing trend from peripheral to central p-Pb collisions. 
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1 Introduction 


Charmonia, bound states of charm and anti-charm quark pairs, are extensively used to study the interplay 
between the perturbative and the non-perturbative regimes of Quantum ChromoDynamics (QCD) iGl] . 
Charmonium production mechanism can be understood as a hard scattering, describing the charm anti¬ 
charm quark pair production, followed by the evolution of the pair into a bound state via a non-perturbative 
process. Models such as colour evaporation (CEM) i,i , colour singlet (CSM) |01 and non-relativistic 
QCD (NRQCD) |Q] are used to describe the charmonium production in hadronic collisions. None of 
these models has so far provided a consistent description of the production cross section and polarisation 
measured in proton-proton (pp) collisions ili. The IS vector state, the J/i/r meson, is abundantly 
produced in hadronic collisions at high energy and measurable through its leptonic decays. Its inclusive 
production contains contributions from direct J/yr, from decays of higher-mass excited states, \I/{2S) and 
Xc, as well as from non-prompt J/t/r, from weak decays of beauty hadrons. 


In proton-nucleus (p-A) collisions, several effects related to the nuclear medium and commonly denoted 
as cold nuclear matter effects (CNM) can affect the production of charmonia. The Parton Distribution 
Eunctions (PDEs) of nucleons bound in nuclei are modified compared to those of free nucleons ffl. 
These functions depend, in particular, on the fraction of the nucleon momentum, Bjorken-x (xbj), carried 
by the probed parton. In the collision energy regime typical of the Earge Hadron Collider (EHC), charm 
quark pairs are produced mainly via the gluon fusion process. The gluon nuclear PDEs (nPDEs) are 
suppressed at low xbj (xbj < 0.01), enhanced at intermediate xbj (0.01 < xbj < 0.3) and suppressed again 
at large xbj (0.35 < xbj < 0.7) compared to those of free nucleons. These three kinematic regions are 
often referred to as the shadowing, anti-shadowing, and EMC regions, respectively. Alternatively, at 
low XBj, the initial colliding nucleus can be described by the Colour Glass Condensate (CGC) effective 
theory Eini] as a coherent and dense (saturated) gluonic system. The kinematic distribution of the 
produced charm quark pairs may be additionally modified by multiple scattering of the incoming gluons 
and/or the quark pairs with the surrounding nuclear medium illl. 11211 or by the energy loss via gluon 
radiation il3l.ll4ll. It was also argued that the interference between the gluons radiated before and after 
the hard scattering can lead to important coherent energy loss effects at large rapidity in the p-going 
direction IlSH . Einally, after their formation, the pre-resonant charm quark pairs or the fully formed 
resonances may interact with the nucleons when passing through the nucleus (nuclear absorption ilhll l 
or with the other particles produced in the p-Pb collision (comovers llTll l. Consequently, they may lose 
energy or fragment into open charm meson pairs. Due to the short time spent by the charm quark pairs 
in the nucleus relative to the J/xf/ formation time at EHC energies, the effect of nuclear absorption is 
expected to be small 


Charmonium production was predicted to be suppressed in a hot medium with a high density of colour 
charges, the Quark-Gluon Plasma (QGP), as a consequence of the colour screening mechanism 0]. 
Such a state can be formed in ultra-relativistic nucleus-nucleus collisions. At the EHC, where the charm- 
quark density is large, charmonium may also be created via the (re)combination of charm quarks either 
during the deconfined phase 12IH or at the phase boundary ll22ll . when the system has cooled down and 
hadronisation takes place. A suppression of J/t/r production in central nucleus-nucleus (A-A) colli¬ 
sions with respect to the one measured in pp collisions scaled by the number of binary nucleon-nucleon 
collisions was observed at the SPS at y^nn ~ 20 GeV |23-25], at RHIC at = 39, 62.4 and 
200 GeV 12^^ and at the EHC at y^^nn = 2.76 TeV 1 33-32]. However, the J/t/r production measure¬ 
ments at the EHC show a much smaller suppression of the yields integrated over transverse momentum, 
Pj, as compared to the results at lower collision energies. The differential results also indicate a smaller 
degree of suppression at px < 3 GeV/c than at higher px ll33l.l34ll. a t mid- and forward rapidity, in agree¬ 
ment with the expectations from (re)combination models j 35 M. Although, qualitatively, the Pb-Pb 
measurements by themselves give a strong indication that the (re)combination effect plays a significant 
role in the J/ t/r production at EHC energies, the quantitative understanding of the involved mechanisms 
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requires a good knowledge of the underlying CNM effects. 


The J/ 1 // production in proton- or deuteron-nucleus collisions was studied at fixed-target (SPS 071.13811 . 
HERA Tevatron and collider experiments (RHIC |41], EHC 1^-3]). At the EHC, a sup¬ 
pression of the J/y production in p-Pb collisions with respect to binary-scaled pp production has been 
observed for pj < 5 GeV/c at large rapidity in the p-going direction and at mid-rapidity, while the 
measurements at high pj as well as at large rapidity in the Pb-going direction are consistent with no 
suppression. The results are in fair agreement with models based on shadowing or coherent energy loss. 
While the J/w suppression at forward rapidity is overestimated by an early CGC calculation 14711 . recent 


calculations Il48l. 14911 are in better agreement with the data. The various CNM effects described above 
should be enhanced at small impact parameters of the collision and thus towards the most central p-Pb 
collisions. Hence, differential measurements as a function of the p-Pb collision centrality are essential 
to further constrain the models, in particular their dependence on the impact parameter of the collision. 

In this paper, we report on new results in p-Pb collisions at = 5.02 TeV for inclusive J/i/r produc¬ 
tion, measured at backward (—4.46 < ycms < —2.96) and forward (2.03 < ycms < 3.53) center-of-mass 
rapidity, ycms> in the decay channel, and at mid-rapidity (—1.37 < ycms < 0.43) in the decay 


channel. Previous measurements have been carried out as a function of rapidity and pj Il42l-l46ll. Here, 
the measurements are performed as a function of the collision centrality, estimated on the basis of the en¬ 
ergy deposited in the Zero Degree Calorimeters (ZDC) At backward and forward rapidity, the J/t/r 
cross section is studied as a function of pj for several centrality classes. The corresponding average val¬ 
ues (pt), and (pj), are extracted from the pj-differential cross sections and the pj broadening, defined 
as A{pj) = (pj)pPb — (Pt)pp’ discussed. The nuclear modification factors are then obtained as a 

function of centrality for the thr'ee rapidity ranges and at backward and forward rapidity, as a function of 
Pj for several classes of centrality. 


2 Detectors and data sets 


The AEICE apparatus and its performance are described in detail in Ref. 115 ill and Ref. Il52ll . respectively. 


Away from the mid-rapidity region, the J/t/r candidates are reconstructed in the decay channel 


using the muon spectrometer ll5lll . covering the pseudorapidity range —4 < T]iab < —2.5 in the laboratory 
frame. The muon spectrometer includes a dipole magnet with an integrated field of 3 T-m, five fracking 
sfafions comprising fwo planes of Cafhode Pad Chambers each, and fwo frigger sfafions consisting of 
fwo planes of Resisfive Plate Chambers each. A sysfem of absorbers is used for filtering ouf fhe hadrons. 
The fronl absorber, made of concrefe, carbon and sfeel wifh a fhickness of 4.1 m (10 nuclear inferacfion 
lengfhs, Aint) is insfalled befween fhe inferacfion region and fhe muon fracking sfafions. A second ab¬ 
sorber, a 1.2 m fhick iron wall (7.2 Aim), is locafed upsfream of fhe frigger sfafions and absorbs secondary 
hadrons escaping from fhe fronf absorber and low-momenfum muons produced predominanfly from n 
and K decays. Einally, a conical absorber placed around fhe beam pipe profecfs fhe specfromefer from 
secondary particles produced in inferacfions of large-T] primary particles wifh fhe beam pipe. 


Af mid-rapidify, fhe J/y candidafes are measured in fhe decay channel wifh fhe cenfral barrel 
defectors in fhe p^udorapidify range |T]iab| < 0.9. The main subsystems used, fhe Time Projecfion 
Chamber (TPC) 15311 and fhe Inner Tracking Sysfem (ITS) i54ll . are placed in a solenoidal magnetic field 
wifh a sfrengfh of 0.5 T. The TPC, fhe main fracking and particle idenfificalion device, is a gaseous drifl 
defector wifh a cylindrical geomefry exfending from 85 fo 247 cm in fhe radial direction and 500 cm 
longifudinally. The particle idenfificalion is performed via fhe measuremenf of fhe specific energy loss, 
dE/dx, in fhe gas volume. 

The ITS, covering a pseudorapidify range |T]iab| < 0.9, consisfs of 6 layers of silicon detectors placed af 
radii ranging from 3.9 fo 43 cm relalive fo fhe beam axis. The fwo innermosf layers are equipped wifh 
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Silicon Pixel Detectors (SPD). The track segments (tracklets) reconstructed from the hits in the two SPD 
layers are used to reconstruct the interaction vertex position and to reject pile-up events (events with 
two or more simultaneous interactions per bunch crossing). The position of the interaction vertex is also 
determined, with better resolution, from the tracks reconstructed in the TPC and the ITS 15ill . 


Two scintillator arrays, VO i55ll . placed on both sides of the interaction point (IP) at —3.7 < T]iab < —1.7 
and 2.8 < T]iab < 5.1, are used as trigger detectors and to remove beam-induced background. They 
are also used for the measurement of luminosity, along with the TO detector 15 ill , consisting of two 
quartz Cherenkov counters, placed on each side of the IP covering the ranges —3.3 < T]iab < —3.0 and 
4.6 < T]iab < 4.9. The Zero Degree Calorimeters (ZDC) located along the beam axis at 112.5 m 
from the IP on both sides, detect protons and neutrons emitted from the nucleus and are used to estimate 
the centrality of the collision. The neutron calorimeter (ZN) is positioned between the two beam pipes 
downstream of the first machine dipole that separates the beams. The proton calorimeter (ZP) is installed 
externally to the outgoing beam pipe. The ZDCs are also used to remove parasitic p-Pb interactions 
displaced from the nominal position. 


The data samples used for the measurements reported in this paper were collected in 2013 in two con¬ 
figurations, obtained by inverting the direction of the p and Pb beams. Due to the asymmetry of the 
energy per nucleon of the p and Pb beams (Ep = 4 TeV and Epb/208 = 1.58 TeV), the nucleon-nucleon 
center-of-mass system is shifted with respect to the laboratory system by Ay = 0.465 in the p-going 
direction. The two beam configurations allow one to measure the J/y/ production in the backward 
(—4.46 < ycms < —2.96) and forward (2.03 < ycms < 3.53) centre-of-mass rapidity (ycms) regions, cor¬ 
responding to the Pb-going and the p-going directions, respectively. They will be further referred to as 
Pb-p and p-Pb, for the first and second case. The dielectron analysis was carried out on a data sample 
corresponding to the p-Pb beam configuration, in the mid-rapidity range —1.37 < ycms < 0.43. 


The dielectron analysis was performed on a sample of events satisfying a Minimum Bias (MB) trigger 
condition and the dimuon analysis used dimuon-triggered events. The MB trigger is defined by a coin¬ 
cidence of the signals from both sides of the VO detector. The efficiency of the MB trigger in selecting 
non-single diffractive p-Pb collisions was estimated to be higher than 99% i57ll . with negligible con¬ 
tamination from diffractive collisions. The dimuon trigger requires, in addition to the MB condition, the 
detection of two unlike-sign muon candidate tracks in the trigger system of the muon spectrometer. This 
trigger selects muons with a transverse momentum px ^ 0.5 GeV/c. This thr'eshold is not sharp in pj 
and the single-muon trigger efficiency reaches a plateau value of ~ 96% at px ~ 1-5 GeV/c. 

In the backward and forward rapidity regions, the measurements are based on a sample of 2.1 x 10^ 
and 9.3 X 10® dimuon-triggered events, respectively. The MB interaction rate reached a maximum of 
200 kHz, corresponding to a maximum pile-up probability of about 3%. The mid-rapidity data sample 
consists of 1.1 X 10^ MB-triggered events, collected at a low interaction rate (~ 10 kHz) and with a frac¬ 
tion of pile-up events lower than 0.6%. Pile-up of collisions from different bunch crossings is negligible 
considering that the 200 ns bunch-crossing spacing is larger than the integration time of the ZDC and 
muon trigger detectors, which are used for the track selection for the dimuon analysis. Two independent 
determinations of the MB trigger cross sections aMB were carried out in the Pb-p and p-Pb configura¬ 
tions using van der Meer scans |58]. The corresponding cross sections amount to = 2.12 ±0.07 b 
and = 2.09±0.07 b, respectively The integrated luminosity is determined as Jf = Nmb/(^mb 

where Vmb is the number of MB events. The number of MB events corresponding to the dimuon- 
triggered sample is evaluated as Vmb = ^ 2 ^/mb ■^DIMU^ where Vdimu is the number of dimuon-triggered 
events and Ti^/mb is the inverse of the probability of having dimuon-triggered events in a MB data 
sample. The determination of Ti^/mb is discussed in Sec. |3l The integrated luminosity was also inde¬ 
pendently measured employing the TO detector. The two luminosity measurements agree within better 
than 1% throughout the whole Pb-p and p-Pb data-taking periods 059l] . The maximum difference is in¬ 
cluded as an additional uncertainty for aMB and thus in the determination of the luminosity uncertainty. 
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The integrated luminosity values used for the results in the backward, forward and mid-rapidity regions 
are 5.81 ±0.20 nb^', 5.01 ±0.19 nb^' and 51.4± 1.9 /rb^\ respectively. 

The centrality determination in Pb-Pb collisions is usually based on charged-particle multiplicity, esti¬ 
mated using the VO signal amplitudes llbOn . However, in p-Pb collisions, the magnitude of the multi¬ 
plicity fluctuations at a given impact pai'ameter is comparable to the whole dynamic range of the MB 
multiplicity distribution. The fluctuations can be related to the various event topologies (e.g. hard colli¬ 
sions with large momentum transfers and/or multiple hard parton-parton interactions, which tend to be 
associated to high multiplicity events, compared to soft collisions without any high-pT particle), detector 
acceptance effects (jets fragmenting in or out the experimental coverage), or other effects, as explained 
in detail in Ref. 0]. Therefore, a centrality selection based on charged-particle multiplicity may select a 
sample of p-Pb collisions that contains biases unrelated to the collision geometry. In contrast, a centrality 
selection based on the energy measured with the ZDC in the Pb-going direction, deposited by nucleons 
produced in the nuclear de-excitation processes following the collision, or knocked out by wounded nu¬ 
cleons, should not induce such biases. The average number of binary nucleon collisions ((A^coii)) or the 
average nuclear overlap function ((TpPb)) for a given centrality class, defined by a selected range of en¬ 
ergy deposited in the Pb-remnant side of ZN, is obtained using the hybrid method described in Ref. 0. 
In this method, the (Acoii) determination relies on the assumption that the charged-particle multiplicity 
measured at mid-rapidity is proportional to the number of participant nucleons ((A^pan))- The values of 
(Apart) for a given ZN-energy class, also noted as ZN class in the following, were calculated by scaling 
the MB value of the number of participant nucleons, (N^), by the ratio of the average charged-particle 
multiplicities measured at mid-rapidity for the considered ZN-energy event classes to the corresponding 
value in MB collisions. The average number of collisions and the average nuclear overlap function were 
then calculated from (Apart) according to the Glauber model ihlll . which is generally used to calculate 


mult\ 


geometrical quantities of nuclear collisions. From here on these values are denoted as and / 

to indicate the ansatz used for their derivation. Other assumptions to derive (Acoii), which are discussed 
in ll50n . use the proportionality of (Acoii) to the yield of high-pj charged particles (10 < pj <20 GeV/c) 
at mid-rapidity or to the charged-particle multiplicity measured with the VO detector in the Pb-going 
direction at forward rapidity. The variations on the (Acoii) values obtained with the three methods do 
not exceed 6% for any of the centrality classes used for this analysis and are taken into account as a 
systematic uncertainty uncorrelated over centrality. Uncertainties of 8% and 3.4% on the determination 
of (A^p) and respectively, are also included as global systematic uncertainties. These uncer¬ 

tainties are obtained by varying the parameters of the Glauber model. Events without a signal in the 
ZN detector, which correspond to very peripheral events 101 . are assigned to the 80-100% centrality 
interval. The values of (A™/^) and (TppjJ^‘) used in this analysi^ are reported in Tab. [H together with 
their uncertainties. 


3 Analysis in the dimuon decay channel 


The analysis approach and the selection criteria are similar to those described in detail in Ref. ll42l] . The 
primary vertex is reconstructed from the hits in the SPD. No specific requirement is applied on the vertex 
properties. In pile-up events, the ZN energy of the two (or more) interactions are summed up, increasing 
the pile-up event contribution in the most central ZN class. This contribution is estimated to be large, 
of the order of 20-30%, for events belonging to the centrality class 0-2% and this class has therefore 
been discarded from the analysis. The muon candidate tracks are reconstructed in the muon spectrometer 
tracking stations using the algorithm described in Ref. 16211 . In order to remove particles at the edge of 
the muon spectrometer acceptance, a fiducial cut on the single-muon pseudorapidity —4 < T]iab < —2.5 
is applied. An additional selection on the radial coordinate of the track at the exit of the front absorber 
(17.6 < Rabs < 89.5 cm) is required to reject muons crossing the high-density section of the absorber. 


*The 0-2% ZN class is excluded in the dimuon analysis due to significant pile-up contribution as detailed in Section[3 
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ZN class 

(Nzn 

/ymult\ 

\^pPb / 

2-10% 

11.7±1.2±0.9 

0.167 ±0.012 ±0.006 

10-20% 

11.0±0.4±0.9 

0.157 ±0.006 ±0.005 

20^0% 

9.6±0.2±0.8 

0.136 ±0.003 ±0.005 

40-60% 

7.1±0.3±0.6 

0.101 ±0.005 ±0.003 

60-80% 

4.3 ±0.3 ±0.3 

0.061 ±0.004 ±0.002 

80-100% 

2.1±0.1±0.2 

0.030 ±0.001 ±0.001 

0-20% 

11.4±0.6±0.9 

0.164 ± 0.009 ±0.006 

60-100% 

3.2 ±0.2 ±0.3 

0.046 ± 0.002 ±0.002 


Table 1: Average numbers of binary nucleon-nucleon collisions and average values of the nuclear overlap 

function (Ip™**) with their uncorrelated and global systematic uncertainty for the used centrality classes. The 
centrality intervals are expressed as percentages of the non-single diffractive p-Pb cross section. 


where energy loss and multiple scattering effects play an important role. Finally, only tracks matching 
the corresponding tracks reconstructed in the trigger stations are selected. 


The J/xj/ candidates are obtained by combining pairs of muons of opposite charge that are reconstructed 
in the rapidity range 2.5 < |yiab| < 4 and with pj < 15 GeV/c. The raw J/t/r yield is estimated for each 
centrality and pj interval from fits of the dimuon invariant-mass distribution performed with various 
functions. For the signal component, an extended Crystall Ball function, which includes non-Gaussian 
tails on either side of the J/xj/ peak, as well as a pseudo-Gaussian function with a mass-dependent 
width 1631] are employed. Due to the poor signal-to-background (S/B) ratio in the tail regions, Monte 
Carlo (MC) simulations are used to constrain the tail parameters in each pj and rapidity interval under 
study. Since there is no degradation of the tracking resolution due to the large occupancy corresponding 
to the most central p-Pb collisions 15211 . the tails are not expected to depend on centrality. The V^(2S) 
resonance is also included in the fit function using the strategy described in Ref. 16411 . For the back¬ 
ground component, two alternative functions are used: a Gaussian with a mass-dependent width and an 
exponential multiplied by a second-order polynomial. The fits are performed in two different invariant 
mass ranges, 2 < m^+^- < 5 GeV/c^ and 2.3 < rn^+^- < 4.7 GeV/c^. In the fitting procedure, the mean 
and width of the J/xj/ signal function, the background parameters and the normalisation factors are left 
free while the tail parameters are fixed fo fhe values esfimafed from fhe simulafions. The obfained J/xf/ 
mass value agrees wifh fhe PDG value 16511 wifhin 5 MeV/c^. The measured widfh increases from 59 
fo 81 MeV/c^ wifh increasing pj. If is found fo be abouf 10% larger fhan in fhe simulafions. The S/B 
ratio in fhe 3a inferval around fhe J/t/r pole increases wifh increasing px and fowards peripheral evenfs, 
ranging from 1.1 (1.3) fo 5.2 (9) in fhe Pb-p (p-Pb) configuration. Figure [J shows examples of fils fo fhe 
unlike-sign dimuon pair invarianl mass dislribulions for fhe Pb-p configuralion for six cenlralily classes 
for px < 15 GeV/c. 


The invarianl mass fils are performed using fhe differenl combinations of signal and background func- 
fions and tiffing ranges described above. The number of J/yr is obfained by inlegrafing fhe J/xj/ signal 
function over fhe filling range. The mean of fhe dislribulion of fhe number of J/tp obfained from fhe 
various fils is used as fhe cenlral value of fhe raw yield, while fhe Roof Mean Square (RMS) is used 
as a sysfemalic uncerlainfy, which ranges befween 0.2% and 3% depending on px and cenlralily class. 
An additional systematic uncerfainly of 2% is added fo fhe signal exlraclion uncerlainfy. If is esfimafed 
from fhe varialion of fhe raw yield when performing fhe til wifh differenl fail paramefers of fhe signal 
funclion. The raw J/xj/ yield varies befween abouf 3000 and 16000 counfs for Pb-p and befween abouf 
5000 and 17000 counfs for p-Pb in fhe cenlralify-differenlial resulfs. In fhe case of fhe cenlralily- and 
Px- double-differential resulfs, if varies befween abouf 70 and 4000 counfs in Pb-p and befween abouf 
150 and 4000 counfs in p-Pb, where fhe lower values correspond fo fhe cenlralily class 80-100% and fhe 
highesl Px range. 
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Fig. 1: (Colour online) Unlike-sign dimuon invariant mass distributions for six centrality classes for pj < 15 
GeV jc in the Pb-p configuration. The solid blue curves correspond to a fit based on a pseudo-Gaussian function 
for the J/y/ and V^(2S) signals (see text) and an exponential multiplied by a second order polynomial function for 
the background. The red dashed lines represent the fitted signal function. 


The J/yr raw yields are corrected for the detector acceptance and efficiency (Axe) estimated from simu¬ 
lations of the J/t/r signal. The muon decay products of the J/are propagated through the experimental 
setup modeled with GEANT 3.21 The procedure used for track reconstruction is the same in data 
and simulations. In the latter, the detector conditions and their variation with time during the data-taking 
period are taken into account. It was checked that the detector occupancy in central collisions does not 
deteriorate the single muon tracking efficiency and resolution, which justifies that only the J/t/r signal is 
simulated and not the underlying p-Pb collision. The pt and rapidity distributions of the J/y signal in 
the simulation were tuned to the reconstructed distributions of the p-Pb and Pb-p data using an iterative 
procedure. The J/y production is assumed to be unpolarised, consistent with the observation that no sig¬ 
nificant J/y polarisation has been measured in pp collisions at = 1 TeV | ^ ^ . The values of A x £ 
integrated over px are 17.1% and 25.4% in the Pb-p and p-Pb configurations, respectively. The lower 
A X £ for the Pb-p configuration is due to a smaller detector efficiency in the corresponding data-taking 
period. The A x £ varies as a function of px from 16% to 33% in Pb-p and from 23% to 48% in p-Pb col¬ 
lisions, where the lowest values correspond to 1 < px < 2 GeV/c and the largest to 8 < px < 15 GeV/c. 
The systematic uncertainty in the choice of the J/t/r kinematic distributions in the simulation is estimated 
by varying the J/ i/r px and rapidity distributions according to the measured ones over various sub-ranges 
of y, Px and centrality (see 1421] for more details). When integrated over px, this uncertainty amounts to 
1.5% for both the Pb-p and p-Pb configurations, while for the px-differential studies it does not exceed 
1.4%. The uncertainty on the dimuon tracking efficiency amounts to 6% (4%) for Pb-p (p-Pb). It is eval¬ 
uated using the difference between the single-muon tracking efficiency obtained from simulations and 
a data-driven approach based on the redundancy of the muon tracking stations 1521] . assuming that the 
efficiencies of the two muons are uncorrelated. This uncertainty is correlated over centrality and is taken 
as constant as a function of px. The uncertainty on the determination of the dimuon trigger efficiency 
has three contributions, which are correlated over centrality. The first uncertainty is due to the statistical 
uncertainty of the trigger detector efficiency which is estimated using data. It is independent of px and 
amounts to 2%. The second uncertainty is extracted from the differences observed between data and 
simulations for the measured trigger response function in the region close to the trigger threshold. This 
uncertainty varies between 0.5% and 3% and is larger at low px- The third uncertainty is due to the small 
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fraction of opposite-sign pairs which are misidentified as like-sign by the trigger system and increases 
from 0.5 to 3% with increasing p^. An additional systematic uncertainty results from the choice of the 
cut, which is applied to the matching of tracks reconstructed in the muon tracking and trigger system. 
Applied to the number of dimuon pairs, this uncertainty amounts to 1% and is correlated over centrality. 

The normalisation factor of dimuon- to MB-triggered events, which is needed to evaluate the 

integrated luminosity, is determined in a two-step procedure as the product where Piii/ip. 

is the inverse of the probability of having dimuon-triggered events (single-muon-triggered 
events) in a corresponding data sample of single-muon-triggered events (MB-triggered events). The var¬ 
ious quantities are estimated from the number of recorded triggered events in each centrality class. This 
factor can also be obtained from the centrality-integrated value scaled by and Adimu/A^dimu’ 

the fraction of MB events and the inverse of the fraction of dimuon events in a given centrality class, 
respectively. The latter method, which is statistically more accurate, is used for the evaluation of F 2 ^/mb- 
The systematic uncertainty is evaluated from the comparison between the two methods and amounts 
to 1 — 2% depending on the centrality class. The value of f 2 /j/MB depends on the centrality class and it 
smoothly increases from 260 and 3340 in Pb-p and from 660 and 3290 in p-Pb from central to peripheral 
collisions. The pile-up event contribution to the number of MB events is estimated by using alternatively 
the interaction vertices reconstructed with the SPD to select pile-up events, or a fast simulation describing 
the ZN energy distribution. The pile-up event contribution is larger in the 2-10% centrality class where 
it amounts to 3.5% and 2.7% in the Pb-p and p-Pb beam configurations, respectively. It decreases to less 
than 2% in all other centrality classes. It has been included in the systematic uncertainty of f 2 /j/MB- It was 
further checked by using the fast simulation that the overall effect of pile-up events, including the shift 
of events from a given centrality class to a more central one, is covered by the systematic uncertainties 
quoted for pile-up events. 


In order to quantify the nuclear effects in p-Pb collisions, reference measurements in pp collisions at 
the same energy are needed. Since there are no experimental data available on the J/y production cross 
section in pp collisions at ^/s = 5.02 TeV, the procedures described in Ref. S for the pt- integrated 
case and in Ref. 14411 for the px-differential case are used. These procedures involve an interpolation 
in energy and an extrapolation in rapidity and are based on existing measurements in pp collisions at 
different energies. The resulting values of the J/i/r cross section interpolated to = 5.02 TeV are also 
reported in those references. 


4 Analysis in the dielectron decay channel 


The analysis method and the selection criteria are similar to those described in detail in Ref. ll44ll . Events 
are selected which contain a primary vertex determined from tracks reconstructed in the TPC and the 
ITS. The vertex position is required to have a distance to the nominal IP smaller than 10 cm along the 
beamline. Due to the lower interaction rate in the data sample used for the mid-rapidity analysis, the pile- 
up event contribution is negligible and the events belonging to the ZN centrality class 0-2% are included 
in the analysis. The electron and positron candidate tracks are reconstructed in the pseudorapidity range 
lOiabI < 0.9. Tracks are required to have at least 70 out of a maximum of 159 clusters in the TPC, a x^ 
normalised to the number of attached TPC clusters smaller than 4 and a distance of closest approach to 
the primary vertex smaller than 3 cm along the beam axis and 1 cm in the plane transverse to the beam 
axis. Only tracks with one associated track point in the innermost layer of the ITS and at least a second 
one in the other layers are selected. This requirement suppresses the secondary electrons from photon 
conversions in the detector material of the ITS. Tracks are required to be compatible within 3.0a with the 
electron hypothesis based on the measured ionisation energy loss of the TPC. In order to reject hadrons, 
the tracks which are compatible within 3.5a with the pion or proton hypotheses are excluded. Tracks 
from identified photon conversions are rejected without any impact on the J/yr signal efficiency. Einally, 
tracks are required to have a transverse momentum larger than 1 GeV/c in order to improve the S/B ratio 
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Fig. 2: (Colour online) Invariant mass distributions of unlike-sign pairs at mid-rapidity for the four con¬ 
sidered centrality classes. The background shape, represented by the dashed blue line, is derived from the event 
mixing technique. The two vertical dashed lines shown in each panel indicate the invariant mass interval used for 
the signal counting. 


in the J/xj/ mass region 1621] . 


The J/xf/ signal is extracted from the invariant mass distribution of candidates. The raw J/t/r 
yields are estimated by bin counting in the invariant mass range 2.92 — 3.16 GeV/c^. The combinatorial 
background is estimated using the event mixing technique, i.e. by pairing electrons and positrons from 
different events. The event mixing is performed in classes of events, sorted according to multiplicity and 
vertex position. The invariant mass distribution of pairs obtained from this procedure is normalised 
to the integral of the same-event unlike-sign dielectron pairs in the mass ranges 2.0 — 2.5 GeV/c^ and 
3.2 — 3.7 GeV/c^, outside of the signal counting interval. A significant fraction of the J/y yield, de¬ 
termined from simulations to be about 30%, falls outside the signal counting window, and is corrected 
for. This is due to the long tail at low masses caused by the bremsstrahlung of electrons in the detector 
material and by the radiative soft photon that may be emitted at the J/y decay vertex. The systematic 
uncertainty on the signal extraction procedure, including uncertainties on the mixed-event background 
scaling and on the invariant-mass shape of the decay channel, is obtained by varying the mass 
region used for the scaling of the mixed-event background and by varying the mass window used for 
counting the signal. Figure [2] illustrates the signal extraction procedure for the centrality classes consid¬ 
ered in this analysis. The central value of the J/y raw yield is obtained as the average of the raw yields 
retrieved from the variation of the signal extraction configurations described above and its systematic 
uncertainty is the RMS of the distribution of the extracted signals. The raw yield has also been evaluated 
using the like-sign method where the residual background after the like-sign subtraction is estimated 
by a linear function and the signal is extracted by bin counting. The two methods have been found to 
provide results compatible within the estimated uncertainties. As a function of centrality, the J/y raw 
yield varies from 73 to 133 counts, with S/B in the interval 2.92 < m^+g- < 3.16 GeV/c^ ranging from 
0.7 to 1.6. 


The Axe correction is estimated with simulations consisting of J/y particles set to decay into an e^e 
pair added to a p-Pb event generated using HIJING Q. The produced particles are subsequently 
propagated thr'ough the experimental setup modeled using GEANT 3.21 and the same reconstruction 
procedure as for the real data is followed. The simulated J/y mesons are assumed to be unpolarised. 
The average value obtained for A x £ is 7.2% with no observed dependence on the collision centrality 
but with a significant dependence on the J/y Pt, having a maximum of ~ 11% at zero pj, a minimum 
of ~ 6% at 2 GeV/c, a second maximum of ~ 10% at 7 GeV/c followed by a slow decrease towards 
higher momenta. This shape is due to the kinematic selections and the momentum dependence of the 
particle identification selection efficiency. 

The uncertainty on the J/y reconstruction efficiency is dominafed by the uncertainty on the electron 
identification. It is estimated from the difference in the TPC specific energy loss distribution of a clean 
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sample of electrons from identified photon conversions in data and electrons from the simulation. For 
dielectron pairs, this uncertainty amounts to 4% and is correlated over centrality. Since the A x £ is not 
constant as a function of pT> its px-integrated value depends on the pj shape used to generate the J/i/r 
mesons in the simulation. To estimate the uncertainty due to the generated px shape, the px-differential 
J/ip spectrum is varied in the simulations. The px distribution is parameterised as 


/(px) =C- 


px 


1 + (pt/po)^)"’ 


( 1 ) 


where C, po and n are parameters constrained by the experimental results in Ref. Il44|] . The systematic 
uncertainty is estimated from the variation of the px-integrated A x £ when varying the fit parameters 
within their uncertainties. It amounts to 3% and is correlated among the centrality classes. Furthermore, 
an uncertainty due to the variation of the spectral shape as a function of centrality is also taken into 
account, given the fact that statistical uncertainties do not allow for a double-differential measurement. 
To evaluate this effect, the relative variation of the fit parameters of the p^ distributions measured at 
forward and backward rapidity for several centrality classes is used. A maximum variation of the A x £ 
by 1.4% is observed and assigned as an uncorrelated uncertainty over centrality. 

The inclusive J/ipproduction cross section in pp collisions at ^/s = 5.02 TeV, which is needed to quantify 
the nuclear effects in p-Pb collisions, is obtained using the interpolation procedure described in Ref. 14411 . 
The method is based on existing measurements in pp collisions at different energies. 


5 Results 

The double-differential J/ ip production cross section for a given centrality class is 


j2^cent 

dydpx 


ycent 

BR 


X O'MB, 


( 2 ) 


where aMB is the p-Pb (Pb-p) MB cross section discussed in Sec. [2l BR is the branching ratio of the 
considered J/ip dileptonic decay channel, which amounts to (5.96 ±0.03)% and (5.97 ±0.03)% for the 
dimuon and the dielectron decay channels 16511 . respectively, and is the inclusive J/ip yield 

per-event. The latter is defined as 


ycent ^ 

. (A X £) • Ay • Apx ’ ^ ’ 

where Aj/y^/+/- is fhe raw number of J/ip mesons decaying info dilepfons for a given cenfralify class, 
rapidify and px range, Vmb is the number of MB events for the given centrality class, A x £ is the 
acceptance times efficiency described in Sec. [3] and |4] and Ay and Apx are fhe widfhs of fhe rapidify and 
Px infervals, respecfively. Table|2]gives a summary of fhe sysfemafic uncerfainfies of fhe J/yr differenfial 
cross secfion, as well as fhe correlations of fhese uncerfainfies over cenfralify, collision sysfem and J/ip 
Px. The px-infegrafed J/ip cross secfions are reported in Tab. [3] for fhe fhree rapidify infervals as a 
funclion of cenfralify expressed in percenfiles of fhe non-single diffractive p-Pb cross secfion. 

Figure |3] shows fhe double-differenfial J/ip cross secfions as a funcfion of px in fhe range 0 < px < 
15 GeV/c af backward (leff panel) and forward (righf panel) rapidify measured for six cenfralify classes. 
The vertical error bars represenf fhe sfafisfical uncerfainfies and fhe open boxes fhe sysfemafic uncer¬ 
fainfies. The sysfemafic uncerfainfies correlafed over cenfralify and px are indicafed as a global relafive 
systematic uncertainly. In order fo characferise fhe evolution of fhe px-differenlial cross secfion shape 
wifh cenfralify, fhe average values (px) and (pj), were exfracfed for each cenfralify class by performing 
a fif fo fhe dala wifh fhe funcfion defined in Eq. [T] The sysfemafic uncerfainfies on fhe dala poinfs fhal 
are correlafed over px are nof considered in fhe til. The uncerfainfies on fhe free paramelers oblained 
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Source of uncertainty 

-4.46 < Jems < -2.96 
cent. (cent, and pj) 

2.03 < Jems ^ 3.53 
cent. (cent, and pj) 

1.37 < Jems 0.43 
cent. 

Signal extraction 

2.0-2.4% (2.8-7.1%) 

2.0-2.1%(2.1-5.3%) 

3.7-7.4% 

p+p^ tracking (I) 

6% 

4% 

- 

p+p- trigger (I) 

3.4% (2.7-3.6%) 

3% (2.7-3.6%) 

- 

p+p^ matching (I) 

1% 

1% 

- 

reconstruction (I) 

- 

- 

4% 

MC input (I) 

1.5% (0.1-1.4%) 

1.5% (0.1-0.4%) 

3% 

MC input 

- 

- 

1.4% 

E2^/MB (III) 

1-3.5% 

1-2.7% 

- 


Uncertainties related to cross section only 


C7mB (I,II,III) 

1.6% 

1.6% 

1.6% 

(Tmb (I,III) 

3% 

3.3% 

3.3% 

BR (I, II, III) 

0.5% 

0.5% 

0.5% 


Uncertainties related to Qpph only 


{T^X) (141,111) 

3.4% 

3.4% 

3.4% 

(7p™“) (II,III) 

1.9-7.2% 

1.9-7.2% 

1.9-5.6% 

^pp (D 

5.3% (8.1-13%) 

5.7% (8.2-11%) 

17% 

CTpp (I, II, III) 

5.5% 

5.5% 

- 


Table 2: Summary of the relative systematic uncertainties for the differential ]/\i/ cross section and Qppb- In the 
backward and forward rapidity intervals, the uncertainties for the px-differential case are indicated in parentheses if 
different from the px-integrated case. Type I stands for uncertainties correlated over centrality. Type II corresponds 
to the uncertainties correlated between the rapidity intervals. Type III is related to forward and backward rapidity 
intervals only and represents the px-correlated uncertainties. The uncertainty on MC input, <7 mb; (Ep™**) and CJpp 
are split into different components according to their correlations over centrality, rapidity intervals and J/xj/ pj. 


ZN class 

da“«/dy(pb) 

-4.46 < yems < -2.96 2.03 < yems < 3.53 

ZN class 

dcT“;^Vdy(pb) 

-1.37 < yems <0.43 

2-10% 

1185±20±49±94 

944±16±33±61 

0-20% 

1582±236±120±98 

10-20% 

1109 ±18 ±32 ±88 

885±14±25±57 

20^0% 

1331±204±55±83 

20^0% 

894±11±23±71 

811±10±21±53 

40-60% 

890 ±160 ±42 ±55 

40-60% 

617±9±15±49 

603±8±15±39 

60-100% 

460±70±52±29 

60-80% 

330±6±10±26 

385±6±11±25 



80-100% 

175±4±4±14 

220±5±5±14 




Table 3: Differential cross sections as a function of centrality. The first quoted uncertainty is statistical while the 
second and third represent the systematic uncertainties, the latter being fully correlated over centrality. 


from the fit of Eq. [Jare propagated to the values of (pj) and (pj). The statistical and systematic un¬ 
certainties on {pj) and {pj) are obtained by performing the fit using, separately, only the statistical or 
the uncorrelated systematic uncertainties on the data points, respectively. The range of integration over 
Pj used to compute {pj) and (pj) is limited to the measured pj interval 0 < px < 15 GeV/c. It was 
verified that extending the integration range to infinity results in an increase of (pj) and (pj) values by 
less than 0.5%. The values of {pj) and {pj) obtained for each centrality class are reported in Tab. H] 
Both (px) and (p|) values increase with centrality, which indicates a hardening of the px distributions 
from peripheral to central collisions in both rapidity intervals. 


In order to quantify the nuclear effects on the Jjxj/pj spectrum shape, the px broadening, A{pj), defined 
as 


MPt) — (PT)pPb “ (Px)pp) 


(4) 


is used. Since there are no measurements for pp collisions at y/s = 5.02 TeV, the value of (Px)pp is eval¬ 


uated from the px-differential cross section in pp collisions calculated with the interpolation procedure 
described in Ref. 11441] . and using the same px-integration range as for p-Pb collisions. Eigure |4] shows 
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Fig. 3: (Colour online) Inclusive J/y/ double-differential cross sections as a function of pj for various centrality 
classes at backward (left) and forward (right) rapidity. The systematic uncertainties correlated over centrality and 
Pi are indicated as a global systematic uncertainty. 


ZN class 

4.46 < Jems 2.96 

{pi) (GeV/c) {p\) (GeV/c2) 

2.03 < 

{pi) (GeV/c) 

Jems 3.53 

{p\) (GeV/c2) 

2-10% 

2.53 ± 0.02 ± 0.03 

9.12 ±0.17 ±0.23 

2.82 ± 0.03 ± 0.03 

11.43 ±0.21 ±0.23 

10-20% 

2.52 ± 0.02 ± 0.03 

9.14 ±0.16 ±0.23 

2.85 ± 0.03 ± 0.03 

11.66 ±0.20 ±0.24 

20^0% 

2.48 ± 0.02 ± 0.04 

8.74 ±0.12 ±0.22 

2.81 ±0.02 ±0.03 

11.36 ±0.15 ±0.23 

40-60% 

2.44 ± 0.02 ± 0.04 

8.48 ±0.14 ±0.22 

2.75 ± 0.02 ± 0.03 

10.98 ±0.16 ±0.24 

60-80% 

2.41 ± 0.03 ± 0.04 

8.41 ±0.20 ±0.22 

2.65 ± 0.02 ± 0.03 

10.22 ±0.18 ±0.21 

80-100% 

2.32 ± 0.03 ± 0.04 

8.03 ± 0.28 ± 0.22 

2.61 ±0.03 ±0.03 

10.00 ± 0.24 ± 0.22 

PP 

2.37 ±0.04 

8.18±0.30 

2.52 ±0.04 

9.28 ±0.40 


Table 4: Values of {pi) and {pj) of inclusive ]/'^f in the range 0 < px < 15 GeV/c. The first quoted uncer¬ 
tainty is statistical while the second is systematic. The values obtained from the pp cross section interpolated to 
= 5.02 TeV are also indicated. 


A(p^) as a function of the number of binary collisions. Our measurements indicate that A(p^) increases 
at backward (forward) rapidity by ~ l.l GeV^/c^ (~ 1.4 GeV^/c^) from peripheral to central p-Pb col¬ 
lisions. At forward rapidity, A(/7j) is larger for all centrality classes and suggests a steeper dependence 
on centrality compared to backward rapidity values. For the most peripheral collisions, corresponding 
to ~ 2, the {pi) value at backward rapidity is compatible with the one in pp collisions, while at 
forward rapidity it is found to be larger than in pp collisions by 0.7 GeV^/c^, which corresponds to 1.4 
times the total uncertainty on the measured difference. The magnitude of the pi broadening observed 
by PHENIX |4l] in d-Au collisions at .^^nn = 200 GeV in the rapidity ranges —2.2 < ycms < —L2, 
bcmsl < 0.35 and 1.2 < ycms < 2.2 is similar to the one measured by AEICE at backward rapidity. At 
forward rapidity, the AEICE data show a stronger pj broadening and a steeper increase with increasing 
centrality as compared to PHENIX results. 


The calculations from Refs. ll72l. l73ll are based on the leading order (EG) GEM production model and 
include initial and final-state multiple scattering of partons with the nuclear medium (denoted as “Mult, 
scattering” in Fig. |4|. The uncertainties on the theoretical calculations are not available. In this model, 
the contribution to pi broadening due to final-state multiple scattering is expected to be sensitive to the 
colour-octet or colour-singlet nature of the pre-resonant cc pair. The calculations are in good agreement 
with the data at backward and forward rapidity. A second model, which is based on a parameterisation 
of the prompt J/y pp cross section and includes coherent energy loss effects from the incoming and 


Ira 


outgoing partons (denoted as “Floss” in Fig. ID), describes well the centrality dependence of A(p|) 
at backward rapidity. The trend predicted by this model at forward rapidity is slightly steeper than 
the data even considering its uncertainty, evaluated by varying the gluon transport coefficient and the 
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< w> 


Fig. 4: (Colour online) pj broadening, A(/7|), as a function of at backward (blue circles) and forward 

(red squares) rapidity. The boxes centered at A(/7y) = 0 represent the total uncertainties of {p\)pp interpolated to 


y/s = 5.02 TeV. The theoretical calculations are from Refs. 07214741] . 


parametrisation of the production cross section. 


In order to study the modification of the J/t/r production in p-Pb collisions with respect to pp interactions, 
the J/xj/ nuclear modification factor is used. For a given centrality class, rapidity and pj range, it is 
defined as 




ycent 

■BR-d’^Opp/dydpT 


(5) 


The (TppjJ*^) values corresponding fo the centrality classes used in this analysis are reported in Tab. [U 
The J/t/r cross section in pp collisions at y/s = 5.02 TeV, d^Opp/dydpj, is obtained by means of the in¬ 
terpolation procedures outlined in Sec.[3]and|4]and described in Refs. |^, The nuclear modification 
factor is usually denoted as RpPb but in this analysis the notation (7pPh is used to emphasise the possible 
bias in the evaluation of (T pjJ*^), as discussed in Sec. |2] and in Ref. iSOll . The systematic uncertainties on 
2pPb are presented in Tab. 


Figure[5]shows the dependence of the px-integrated 2pPb on the collision centrality, expressed as (A^^p')- 
The results for the backward, mid- and forward rapidity intervals are displayed in the left, middle and 
right panel, respectively. At backward rapidity, the 2pPb values show that the measured J/ tp production 
is compatible, within the total uncertainties, with expectations from binary collision scaling for all cen¬ 
trality classes. When considering only the uncertainties that are not correlated over centrality, an increase 
from peripheral to central p-Pb collisions is observed in the data. At forward rapidity, the J/yr yield is 
suppressed with respect to the binary-scaled pp reference for all the considered centrality classes. The 
values of 2pPb measured at forward rapidity exhibit a decrease from 0.85 for the 80-100% centrality 
class down to 0.66 for the 2-10% centrality class. Within the present uncertainties, the mid-rapidity 
results suggest a similar degree of suppression of the J/yr yield as at forward rapidity and no conclusion 
can be drawn on a possible centrality dependence. 


Our measurements are compared to several theoretical models including a next-to-leading order (NFO) 
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Fig. 5: (Colour online) Inclusive J/y/ gppb as a function of at backward (left), mid (center) and forward 

(right) rapidity. The boxes centered at Qpp b = 1 represent the relative uncertainties correlated over centrality. The 
theoretical calculations are from Refs. Ill6lll7ll741751] . 


CEM calculation |[l^. 75] which contains the EPS09 NEO nPDE paramctcrisation [0] (denoted as 
“CEM+EPS09 NEO” in Eig. |5]l, a model employing the EPS09 EO nPDE with or without effects from 


the interaction with a comoving medium llTI] (denoted as “EPS09 EO+comovers” in Eig. O, and the 
coherent energy loss model 17411 described above. In the CEM+EPS09 NEO and EPS09 EO+comovers 
models, assuming the J/xf/ production process is gg —)■ J/t/r (2 —1), the ^bj values of the gluon from the 
Pb nucleus span a range of about 1 • 10^^ < vbj < 5 • 10^^ at backwai'd rapidity, 4 • 10^"^ < xbj < 2 • 10^^ 
at mid-rapidity, and 2 • 10^^ < xbj < 8 • 10^^ at forward rapidity. The backward rapidity interval there¬ 
fore corresponds to the xbj range in the transition between the anti-shadowing and the shadowing re¬ 
gion, whereas the mid- and forward rapidity intervals probe a region for which the gluon shadowing is 
expected to be strong. The CEM+EPS09 NEO model uncertainties are evaluated from the EPS09 un¬ 
certainty, which gives the dominant contribution, and from a variation of the values of the charm quark 
mass, the normalisation and the factorisation scales in the pQCD calculation. The CEM+EPS09 NEO 
model reproduces well the centrality dependence in each rapidity range. At mid- and forward rapidity, 
the data are better reproduced when a strong shadowing is considered in the model. In the framework 
of the EPS09 EO+comovers model, the presence of a comoving medium has only a small effect onJ/xj/ 
production at forward rapidity since its density is expected to decrease towards the p-going direction. 
The effect of comovers is more pronounced at mid-rapidity and especially at backward rapidity and it 
increases with increasing centrality. The uncertainties on these theoretical calculations are not available. 
At backward rapidity, the increase of 2pPb towards central collisions observed in the data is better re¬ 
produced when the comover effect is not included in the model. Einally, the shape and magnitude of 
2pPb is well described by the Eloss model in all rapidity intervals, although the model does not predict 
an increase with increasing centrality at backward rapidity, as indicated by the data. 


surements also include the contribution of J/xj/ mesons from Zj-hadron decays. The 
tracted from using the relation 2^™ 
non-prompt to prompt J/xj/ production cross sections and 2pPb 


It is worth pointing out that the calculations above are done for prompt J/xj/ production, while the mea- 

can be ex- 

= QpPb+fs ■ (Gppi - Gppr'’™'"'’')' where fe is the ratio of 

is the nuclear modification factor 
of non-prompt J/xj/ mesons. A value of /b of about 0.11 at 2 < yems < 4.5 and for pj < 14 GeV/c 
can be calculated from the EHCb measurements in pp collisions at ^/s = 1 TeV 17611 . The value of /g 
does not show a strong variation within the quoted rapidity range and with energy, as indicated by the 
comparison with the results in pp collisions at ^/s = 8 TeV i77n . Hence, the value of /g calculated 
at -y/? = 7 TeV in 2 < yems < 4.5 is used for the following. At mid-rapidity, a value of fs of about 
0.17 at |y| < 0.9 and integrated over pj can be extracted from the measurements of AEICE in pp col¬ 
lisions at = 1 TeV 17811 . The nuclear modification factor of non-prompt J/t/r was measured to be 
0.98 ± 0.06 ± 0.01 (0.83 ± 0.02 ± O.O^for -4 < y^ms < -2.5 (2.5 < yems < 4) and px < 14 GeV/c at 
y/^NN = 5.02 TeV in p-Pb collisions i43n . If the non-prompt f/xp 2pPb> which has not been measured as 
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Fig. 6: Inclusive J/t/ Qppb as a function of pj for the 2-10%, 10-20%, 20^0%, 40-60%, 60-80% and 80-100% 
(from top to bottom) ZN centrality classes at backward (left) and forward (right) rapidity. The boxes centered at 
CpPb = 1 represent the relative uncertainties correlated over pj. The theoretical calculations are from Refs. 0 
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a function of the centrality, is conservatively assumed to vary from 0.6 to 1.3 in each centrality interval, 
then the differences between the inclusive and prompt J/yr nuclear modification factors cannot exceed 
15% in any of the centrality classes and are smaller than the quoted uncertainties. 


Figure 0 presents the pj-dependence of for the 2-10%, 10-20%, 20-40%, 40-60%, 60-80% and 
80-100% centrality classes, from the top to the bottom panels, respectively. The left (right) panels show 
the backward (forward) rapidity results. At backward rapidity and for the most central collisions (2-10% 
and 10-20% centrality classes), 2pPb is compatible with unity in the full pT interval, and an increase of 
2pPb from pt < 1 GeV/c to pj > 1 GeV/c is suggested by the data. For semi-central and peripheral 
collisions, 2pPb is compatible with unity over the full pT range. At forward rapidity, for all centrality 
classes other than the most peripheral one, J/ ip production is suppressed compared to binary-scaled pp 
production at low px. For these centrality classes, 2pPb increases with increasing px and is compatible 
with unity at high px within the uncertainties. The magnitude of 2pPb slightly increases from central to 
semi-peripheral collisions over the full px range. For the most peripheral collisions 2pPb is compatible 
with unity and, with the current uncertainties, it does not show a significant px dependence. 


The data are compared to the calculations from the CEM-I-EPS09 NFO S F/sh and the Eloss iF^ 
model. The CEM-I-EPS09 NFO calculations describe reasonably well the 2pPb results at backward and 
forward rapidity. The Eloss model reproduces well the p^ dependence of 2pPb at backward rapidity 
for ah centrality classes. At forward rapidity, a good agreement of the calculations with the data is 
observed for peripheral collisions (60-80% and 80-100% centrality classes), while the _ 2 j dependence 
becomes steeper than in data towards more central collisions. It was observed in Ref. 14411 that a better 
agreement is reached with the px dependence of the nuclear modification factor in centrality-integrated 
p-Pb collisions at forward rapidity when shadowing effects are included in the model. 


6 Conclusions 


The cross sections and nuclear modihcation factors, 2ppb, of inclusive J/ip production have been mea¬ 
sured with AFICE as a function of rapidity, px and centrality in p-Pb collisions at .^^'nn = 5.02 TeV. 
For the most peripheral p-Pb collisions, no modification with respect to pp collisions is observed within 
the uncertainties of the measurements for both the shape of the px spectrum of the J/y and the 2pPb 
measurements. On the contrary, the results in central p-Pb collisions suggest sizeable nuclear effects. At 
both backward (Pb-going direction) and forward (p-going direction) rapidity the A{pj) measurements 
show a Px broadening which increases monotonically from peripheral to central p-Pb collisions with 
larger values at forward rapidity. Our measurements show a stronger px broadening and a steeper in¬ 
crease with increasing centrality at forward rapidity as compared to PHENIX results in d-Au collisions 
at = 200 GeV 14111 . At backward rapidity, a modest increase of J/yr production compared to a 
binary-scaled pp reference is suggested by the data in most central collisions. At mid-rapidity, the data 
indicate that J/ ip production is suppressed compared to binary-scaled pp cross section over the entire 
centrality range. Within the current uncertainties, the increasing suppression towards central p-Pb colli¬ 
sions suggested by models is compatible with the data. Finally, at forward rapidity, a clear suppression, 
which increases towards central events, is observed. The px- and centrality-differential results show that 
the suppression is stronger at low px and tends to vanish at high px. Given the uncertainties of both 
the measurements and the theoretical calculations, we observe a fair agreement of the models based on 
coherent energy loss and multiple scattering with the measured px broadening. Models based on nPDF 
and coherent energy loss are in fair agreement with the nuclear modification factor measurements. The 
results presented in this paper provide an important baseline for understanding and constraining the cold 
nuclear matter effects in p-Pb collisions as well as their centrality dependence. Such an information is 
essential for a quantitative interpretation of the results obtained in Pb-Pb collisions. 
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